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We selected 40 candidate Lyman Alpha Emitting galaxies (LAEs) at z — 3.1 with observed frame equivalent 
widths > 150A and inferred emission line fluxes > 2.5 x 10~ 17 ergs cm" 2 s" 1 from deep narrow-band and broad- 
band MUSYC images of the Extended Chandra Deep Field South. Covering 992 arcmin 2 , this is the largest 
"blank field" surveyed for LAEs at z ~ 3, allowing an improved estimate of the space density of this population 
of 3 ± 1 x 10~ 4 /i7 Mpc~ 3 . Spectroscopic follow-up of 23 candidates yielded 18 redshifts, all at z — 3.1. Over 
80% of the LAEs are dimmer in continuum magnitude than the typical Lyman break galaxy spectroscopic limit 
of R = 25.5 (AB), with a median continuum magnitude R ~ 27 and very blue continuum colors, (V -z) — 0. 
Over 80% of the LAEs have the right UVR colors to be selected as Lyman break galaxies, but only 10% also 
have R < 25.5. Stacking the UBVRIzJK fluxes reveals that LAEs have stellar masses ~ 5 x lO 8 ^ 2 , M and 
minimal dust extinction, Ay < 0.1. Inferred star formation rates are ~bh~$s/\.Qyx~ l , yielding a cosmic star 
formation rate density of 2 x 10~ 3 /i7oMQyr _1 Mpc~ 3 . None of our LAE candidates show evidence for rest- 
frame emission line equivalent widths EW„,„ >240A which might imply a non-standard IMF. One candidate is 
detected by Chandra, implying an AGN fraction of 2 ± 2% for LAE candidate samples. In summary, LAEs at 
Z ~ 3 have rapid star formation, low stellar mass, little dust obscuration and no evidence for a substantial AGN 
component. 1 

Subject headings: galaxies:high-redshift 



1. INTRODUCTION 

Because Lyman a emission is easily quenched by dust, the 
Lyman Alpha Emitting galaxies (LAEs) are often character- 
ized as protogalaxies expe riencing their first burst of star for- 
mation l lHu & McMahonll 19961) . However, the differing be- 
havior of Lyman a and continuum photons encountering dust 
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and neutral gas makes it possible for older galaxies to exhibit 
Lyman a emission when morphology and kine matics favo r 
the escape of these photons (e.g. Haiman & Spaans 1999). 
Hence the LAEs could instead represent an older population 
with actively star-forming regions. 

LAEs offer the chance to probe the bulk of the high-redshift 
galaxy luminosity function as the strong emission line allows 
spectroscopic confirmation of objects dimmer than the contin- 
uum limit R < 25.5. Previous studies of LAEs at z ~ 3 have 
concentrated on known overdensities ij Ste idel et alJ 120001 
Havashino et al. 20041 IVenemans et al.ll20Q5l> or searches for 
Lyman Alpha emiss ion near known D am ped Lyman Alpha 
absorption systems ( Fvn bo et aLll2003l see IWolfe etaLll2005l 
for a review). Blank fields, i.e. those not previously known 
to contain unusual objects or overdensities, ha ve been s tudied 
at z = 3.1 and z = 3.4, covering 468 arcmin 2 jC iardullo et alJ 
120021) and 70 arcmin 2 (ICowie & Hull998tlHu et al.ll998l) . re- 
spectively. Significant work has been done in recent years 
on large blank fields at higher redshi fts to study the LAE 
luminosity func tion at z = 3.7 (| FujrtaeLffl I2QQ3I ), z = 4.5 



JHu et alJll998T) z = 4.9 (lOuchi et al.H2003l I Shimasak u et al 
1 20031) . z = 5.7 (Martin & Sawickil l2004lMalhotra & Rhoads 
2004 and references therein) and z = 6.5 (Malhotra & Rhoads 

2004) . Spectroscopically con firmed samples are sm all, in- 
cluding 31 LA Es at z = 3.1 (IVenemans et al J 120051). 18 at 
z = 4. 5 (LALA, iDawson et al.N2004l) 27 at z = 5.7 (iHu et alJ 
| 2004[ lOuchietal] 120051) . and 9 at z = 6.6 ( iTaniguchi et alJ 

2005) . T he current invest i gatio n expands upon the blank-field 
survey of lCiardullo et alJ ( 120021) by covering twice the area to 
a narrow-band detection limit one magnitude deeper. 

The study of Lyman Alpha Emitting galaxies at 
z ~ 3.1 is a major goal of the Multi waveleng t h Sur- 
vey by Yale-Chile (MUSYC, Gawiser eta!] 120061 
|http : / /www . astro . yale . edu/MUSYC) . The Ex- 
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tended Chandra Deep Field South (ECDF-S) has been 
targeted with deep narrow-band imaging and multi-object 
spectroscopy, complemented by deep broad-band UBVRIzJK 
and public Chandra+ACIS -I imaging. These multi wavelength 
data make it possible to study the physical nature of LAEs and 
to distinguish star formation from AGN as the source of their 
emission. We assume a ACDM co smology consistent with 
WMAP results dBennett et alJl2003l) with Q, m = 0.3, ft A = 0.7 
and Hq = IQh-iQ k m s" 1 Mpc" 1 . All ma gnitudes are given in 
the AB95 system (Fukug itaet al.ll9 96). 

2. OBSERVATIONS 

Our narrow-band imaging of the ECDF-S was ob- 
tained using the NB5000A filter (50A FWHM) with 
CTI04m+MOSAIC-II on several nights from 2002 to 
2004 for a total of 29 hours of exposure time. Our 
UBVRI imaging results from combining public images 
taken with ESQ2.2m+ WFI by the ESQ Deep Public Survey 
and COMBO- 17 tea ms ( Erbe n et alJl2005l iHildebrandt et al] 
120051 lArnouts et aljl2001t IWolf et alJl2004l) . Our 7 imag- 
ing was taken with CTI04m+MOSAIC-II on January 15, 
2005. Details of our optical images will be presented in 
E. Gawiser et al. (in prep.). Our JK images of ECDF-S 
were obtained with CTI04m+ISPI on several nights during 
2003-2004 and will be described in E.N. Taylor et al. (in 
prep.). The final images cover 31.5' x 31.5' = 992 arcmin 2 
centered on the Chandra Deep Field South and were pro- 
cessed through the MUSYC photometric pipeline to gener- 
ate APCOR R (corrected aperture) fluxes and uncertainties as 
described in iGawiser et al J (12006). Table[2gives our source 
detection depths. 

TABLE 1 

5cr Point Source Detection Limits for MUSYC ECDF-S Images 
in AB magnitudes. 



NB5000 U B V R I z' J K 
25.5 26.0 26.9 26.4 26.4 24.6 23.6 22.7 22.0 



Multi-object spectroscopy of 23 LAE candidates was per- 
formed with the IMACS instrument on the Magellan-Baade 
telescope on Oct. 26-27, 2003, Oct. 7-8, 2004, and Feb. 4-7, 
2005. The 300 line/mm grism was used with 1.2" slitlets to 
cover 4000-9000A at a resolution of 7.8A . Details of our 
spectroscopy will be given in P. Lira et al. (in prep). 

3. CANDIDATE SELECTION 

The greatest challenge in selecting Lyman Alpha Emitting 
galaxies at z — 3 . 1 is to minimize contamination from z — 0.34 
galaxies exhibiting emission lines in [O II] 3727 A . These in- 
terlopers can be avoided by requiring a high equivalent width 
(>150A in the obser ved frame) which el iminate s all but th e 
rarest [O II] emitters JTerlevich et alJl991tlSte"rn et alJ2000l) . 
Contamination from [O III]5007A is minimal at these wave- 
lengths, as the volume for extragalactic emitters is small, and 
Galactic planetary nebulae are very rare at such high Galactic 
latitude (b = -54). 

Selecting LAEs requires an estimate of the continuum at the 
wavelength of the narrow-band filter, so we tested weighted 
sums of the B and V flux densities and found that f^ v = 
(fS +f„)/2 minimizes the scatter in predicting the NB5000 
flux density of typical objects. The "narrow-band excess" 
in magnitudes, BV - A/65000, was then used to select the 



candidate LAEs with (BV -NB5000) > 1.5, corresponding 
to EW^ a > 150A . When the broad-band fluxes are small, 
significant errors in the equivalent width estimate may result, 
and a small fraction of the numerous objects without emission 
lines, i.e. with (BV -NB5000) ~ 0, could enter the "narrow- 
band excess" sample. To avoid both types of interlopers, we 
calculated a formal uncertainty in the BV -NB5QQ0 color in 
magnitudes, (Jbv-nb, and required candidate LAEs to have 
(BV-NB5000)-a BV - NB > 1.5 and (BV-NB5000)-3a BV . NB > 
0. The latter crit erion i s similar to the color excess re- 
quirement of B unker et alJ ill 9951) . but our color uncertainties 
are object-specific a nd account for varia tion in image depth 
across the field (see IGawiser et al]l2006l) . To make spectro- 
scopic confirmation feasible, we also required NB5000<25.0, 
implying an emission line flux > 2.5 x 10~ 17 ergs cm -2 s _1 . 
Visual inspection to eliminate false narrow-band detections 
caused by CCD defects or cosmic ray residuals resulted in 
40 candidate LAEs. 23 of these candidates have been ob- 
served spectroscopically, yielding 1 8 confirmations where the 
Lyman a emission line was clearly detected in both the two- 
dimensional and extracted spectra and no other emission lines 
were visible across the full optical spectrum. We tested the 
procedure by observing lower-equivalent width objects and 
found several [O II]3727 emitters; all of these interlopers ex- 
hibit clear emission lines in H/3, [O III]4959,5007 and Ha. 
Five of the LAE candidate spectra fail to show emission lines. 
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FIG. 1. — Histogram of R-band magnitudes for candidate LAEs (thin 
histogram) and the subset of confirmed LAEs (thick histogram), with typical 
spectroscopic limit of R = 25.5 marked with dashed vertical line. Objects 
with negative R fluxes were assigned R = 30. 

4. RESULTS 

Figure ^ shows the distribution of candidate and con- 
firmed LAE 7?-band continuum magnitudes versus the "spec- 
troscopic" Ly man break galaxy (LBG) limit of R < 25.5 
JSteidel et al.ll2003l) . Our study of LAEs is able to observe 
objects much dimmer than this, with a median magnitude 
R ~ 27. 36/40 candidates and 15/18 confirmed LAEs have 
R > 25.5, showing the efficacy of LAE selection in identi- 
fying objects from the bulk of the high-redshift galaxy lu- 
minosity function. Figure |2] shows UV cor rR colors of our 
LAE candidates v ersus t he LBG selection region determined 
by IGawiser et alJ (|2006), where V corr refers to the V-band 
magnitude after subtracting the flux contributed by the Ly- 
man a emission lines. Only 2 out of 18 confirmed LAEs 
fall within the R < 25.5 "spectroscopic" LBG sample, but 
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16 out of 18 fall within the color selection region. About 
half of our candidate LAEs would meet the R < 27 mag- 
nitude limit of the "photom etric" LBG sample explored by 
ISawicki & Thompsonl (12005), and these objects should com- 
prise 5% of their sample. 
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FIG. 2. — UVR color-color plot of confirmed LAEs (solid circles), can- 
didate LAEs with spectroscopy but no confirmed redshift (open circles) and 
candidate LAEs without spectroscopy (plusses) versus distribution of the en- 
tire 84,410 object optical catalog (dots). The polygonal region in the upper 
left is the Lyman break galaxy selection region. 

In order to investigate the full SED of the LAEs, which 
are too dim to obtain individual detections in our NIR 
photometry, we measured stacked fluxes for the confirmed 
sa mple and show the result s of SED modelling in Fig. 
|3] iBruzual & Chariot! (|2003) population synthesi s models 
were used, with constant star formation rate, a |$alpeter 
( 119551) initial mas s function from O.l Mm to 100M Q , so- 
lar metallicity and lCalzetti et alj (|1 997|) dust red dening (e.g. 
iForster Schreiber et alJl2004t Ivan Dokkum et alj r2004). Un- 
certainties in the stacked photometry were determined using 
bootstrap resampling and are close to the formal errors cal- 
culated from the reported APCORR flux uncertainties. Pa- 
rameter uncertainties were computed via a Monte Carlo anal- 
ysis where the stacked fluxes were varied within their un- 
certainties to yield a probability distribution of best-fit pa- 
rameters. The age of the stellar population is weakly con- 
strained and has been restricted to the physically reasonable 
range 10 Myr < < 2 Gyr. The best-fit parameters shown 
in Fig. |3] correspond to minimal dust extinction, significant 
star formation rates (5<SFR<23 h% M Q yr _1 at 95% confi- 
dence) and low stellar mass (the 95% confidence upper limit 
is M* = 8.5 x 10 9 /z7qMq). The LAEs appear to have much 
less dust and stellar mass than the ~ 500 My r old, Ay ~ 1, 
~2x 1O 1O M0 Lyman break galaxy population ( Shaplev et al. 
120011) or the ~ 2 G yr old, A v ~ 2.5, ~ lO^M^ D istant Red 
Galaxy population (Forster Schrei ber et alj [2004). The star 
formation rates of the confirmed LAEs inferred from their 
Lyman a luminosities average 5hjl M Q yr _1 and from their 
rest-frame UV continuum luminosity densities average 9hjf } 
M Q yr _1 . The consistency of these values with the best-fit SFR 



from SED modelling implies minimal dust extinction. 
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FIG. 3. — UBVRIzJK broad-band photometry (average flux density of 
stacked sample) of confirmed LAEs along with best-fit model from SED fit- 
ting (solid) with model parameters listed. The dotted curve shows a maxi- 
mally old model with stellar population age fixed to 2 Gyr (the age of the 
universe at z = 3.1), A v = 0.1, SFR=lh^ M yr 1 and M* = 1.1 X lO 10 /^ 
Mq. 

To check for AGN contamination of our LAE candidate 
sample, we have looked for Chandra detections of these ob- 
jects. O ne LAE candidate h as a n X-ray detect i on in t he cat- 
alogs of lVirani etail ll2005l) and|Eehmer et al. (2005af), with 
a 0.5-8keV luminosity of 10 44 erg s -1 . No other candi- 
dates showed individual detections, so we removed this ob- 
ject and performed a stacking analysis (e.g. lRubin et all2004t 
lLehmeretall2 005b) which resulted in a non-detection of the 
entire population. Using the conversio n between SFR and 
X-ray flux given by Ranalli et al. (2003), the upper limit on 
the average star formation rate per object is 200^1 M Q yr" 1 , 
which is clearly consistent with the observed SFR. None of 
our LAE spectra show broad emission line widths (> 1000 
km s" 1 ) that would be inconsistent with the energetics of star 
formation. We therefore expect that very few LAE candi- 
dates contain luminous AGN which dominate their Lyman a 
or continuum emission. 

5. DISCUSSION 

Our survey covers 31.5' x 31.5' x (Az = 0.04) or 59 x 
59 x 38/i7QMpc 3 , yielding an LAE number density of 3 ± 1 x 
10~ 4 /z 3 Mpc~ 3 , equivalent to 4000 ± 1600 deg" 2 per unit red- 
shift. The survey volume was computed using the filter band- 
pass FWHM=50A, and the five candidates without confirmed 
redshifts were assumed to be LAEs. The error bars account 
for variations in the LAE abundance within our survey vol- 
ume caused by large-scale structure assuming a bias of 2. The 
true uncertainties could be bigger given the large fluctuations 
in den sity observed for LAEs at z = 4.9 bv fShimasak u et all 
(2004). Combining the measured number density and using 
the best-fit star formation rate per object of 6/1^0 M Q yr" 1 , we 
find a cosmic star formation rate density of 2 x 10~ 3 /i7o M Q 
yr _1 Mpc" 3 . This is significantly less than the LBG SFR den- 
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sity JSteidel et al.lll999l) . but it underestimates the total LAE 
contribution due to our requirements of high equivalent width 
and relatively bright AfB5000 flux designed to select a pure 
sample amenable to spectroscopic confirmation. A detailed 
calculation of the LAE luminosity function at z — 3.1, which 
can be integrated to give a fuller estimate of the SFR density, 
will be given in C. Gronwall et al. (in prep). 

The number density, stellar masses, star formation rates, 
and median U V continuum fluxes fo und for LAEs are within 
a fact or of three of those predicted bv lLe Delliou et aD J2005t 
2006); the agreement is even better when our equivalent width 
threshold is accounted for. The only strong disagreement seen 
versus these models is their claimed escape fraction of 0.02 
for Lyman a photons versus our lower limit of 0.2 (and best- 
fit of 0.8) implied by the comparison of star formation rates 
determined from the observed Lyman a luminosities and SED 
modelling. This discrepancy could be resolved by using a 
larger escape fraction and a standard IMF instead of the top- 
heavy IMF assumed in the models. 

Our determination that z = 3 . 1 LAEs are pre dominantly 
blue contrasts with th e results of Stiavelli et al. (200jJ and 
iPascarelle et all (119981) that LAEs in blank fields at z ~ 2.4 
are typically red, (B-/)~ 1.8. This differs from the me- 
dian value of (V con -z) — 0.1 for our spectroscopically con- 
firme d LAEs and the median color (V-I) ~ 0.1 measured by 
IVenemans et alJ fc005l) . The difference seems unlikely to be 
caused by evolution in the LAE population from z = 3.1 to 
z = 2.4 given the sma ll increase in the age of the universe. 

At z = 4.5, LALA JMalhotra & R hoads 2002) reported that 



a majority of LAE candidates had EW, ej , > 240A , providing 
evidence of a top-heavy IMF possibly caused by Population 
III stars, although equivalent widths this high could also re- 
sult from highly anisotropic radiative transfer due to the dif- 
fering effects of dust and gas on Lyman a and UV-continuum 
photons. This photometric measurement is sensitive to con- 
siderable scatter when the sample is selected in the narrow- 
band and the broad-band imaging is shallow, as broad-band 
non-detections can receive extremely large implied equiva- 
lent widths, and this is guaranteed to occur for any spurious 
narrow-band detections. Indeed, when 2a upper limits on the 
continuum flux were used, only 10% of their z = 4.5 sample 
had such high EWs, and ~ 20% of the confirmed objects have 
EW res , > 240 A jDawson et all2004l) . We do not find equiva- 
lent widths this high for any of our LAE candidates at z = 3.1 . 
The difference might reveal evolution in the LAE population 
or could be the result of small number statistics. 
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